ANALYTIC STUDY OF PERSISTENT CURRENT IN A 
TWO-CHANNEL DISORDERED MESOSCOPIC RING 

Jean HeinrichJ*] 

Departement de Physique, B5a, Universite de Liege, 
Sart Tilman, B-^OOO Liege, Belgium 
(Dated: January 11, 2013) 

Abstract 

We present an extensive analytical study of persistent current in a weakly disordered two-chain 
cylindrical ring threaded by an Aharonov-Bohm flux < <p < <po/2 (with 4>q the flux quantum) 
and described by the Anderson model. The effect of the disorder reveals a strong reduction of the 
persistent current for flux values near </>o/4. 

In conjunction with the pure system (zeroth order) current profile averaged over numbers of 
electrons and earlier results for the effect of disorder in one-dimensional rings, our two-channel 
results provide a simple interpretation of salient features of numerical results of Bouchiat and 
Montambaux (BM) for persistent current in an assembly of many-channel disordered rings. Single- 
channel (one-dimensional) effects are responsible for the dip in the persistent obtained by BM near 
= and the corresponding peak near (po/2, while the effect of disorder in independent channel 
pairs accounts for abrupt decreases of current superimposed to a continuous linear decay as the 
flux value 0o/4 is approached from above and from below, respectively. The persistent current in 
the two-channel ring involves a free particle current averaged over electron numbers of periodicity 
4>0 /2, and a dominant disorder effect which has periodicity 4>q. 
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I. INTRODUCTION 



An Aharonov-Bohm flux threading a metallic or semiconducting ring leads to a persistent 
equilibrium current even if the ring is disordered, provided it is of mesoscopic size smaller 
than the mean free path for inelastic scattering jj. This is a consequence of the quantum 
coherence of the electronic wavefunction, which is not destroyed by elastic impurity scatter- 
ing. Intensive study of persistent current was initiated many years ago with the appearance 
of the seminal paper by Biittiker, Imry and Landauer^ predicting persistent current in a 
one-dimensional disordered ring. The persistent current is periodic in the magnetic flux <ft 
with a period expressed in terms of the flux quantum O — hc/e (with h the Planck constant, 
c the speed of light and — e the electron charge). The periodicity of the persistent current is 
due to the modification of the periodic boundary condition for the wavefunction ip{x) along 
the ring of length L in the presence of the flux 0, namely 

vp(x + L) = ip(x)e . (1) 

The persistent current in a particular energy level E n is related to the corresponding flux 
derivative by 



J « = ~ c ^nr ' ( 2a ) 



where c is the speed of light. The total persistent current in the ring is 



(2b) 



The access to quantum-mechanical features of solids renders accurate measurements of 
persistent current highly desirable. An important step in this direction has been achieved 
recently with the a_ppearance of high-precision measurements of persistent current by two 
different groups {jj, {J showing for the first time quantitative agreement with theoretical 
models for non-interacting electrons in diffusive rings (see also p, |6]). 

The work of Biittiker et al. 2\ for a one-dimensional ring was later amplified and extended 
by detailed analytical studies p], |8| of the effect of disorder on persistent current, using the 



Anderson tight-binding model with random atomic site energies. We focus, in particular, 
on the analysis of js| (referred to as I in the following) which leads to a convenient exact 
perturbation theory for the energy levels of the ring, which will be generalized below in 
the case of a disordered two-channel ring. Persistent current has been discussed in [7, s| 
in various flux intervals close to <fi = 0. Throughout this paper we consider specifically the 
interval < < <po/2 in which extensive numerical simulations of persistent current have 
been performed by Bouchiat and Montambaux (BM)[9] for ten-channel disordered rings. 
An important objective of our work is to interpret the simulation results[9| by adapting the 
results of for a one-dimensional (one-channel) ring and new results for persistent current 
in a two-channel disordered ring (the thinnest multi-channel ring!) whose derivation forms 
the core of the present article. The ring-model studied numerically in [9] is a quasi-one- 
dimensional Anderson model with on-site disorder and nearest-neighbour inter- as well as 
intra-chain hopping. We focus attention on the fig. 4 of 9] which reveals, in particular, 
the existence of an important dip in the persistent current near = and a corresponding 
peak near (fio/2. Abrupt current decays are also observed in the figure when approaching 
the value 0o/4 from below and from above, respectively. Another specific feature of the fig. 
4 in [9j is that the persistent current averaged over even- and odd numbers of electrons 10 ] 
has periodicity 0o/2. An important observation of BM is that the shape of their fig. 4 is 
independent of the number of channels, M. This suggests indeed that important features of 
the figure might be observed already in single- and two-channel systems as is confirmed in 
Sect llVBl below, by adapting, in particular the results of I (averaging over even- and odd 
numbers of electrons). 

The special interest in studying transport in disordered quasi-one-dimensional mesoscopic 
rings involving a large number M, of parallel channels, is the existence in these systems of 
a wide domain of metallic conduction extending over length scales between the elastic mean 
free path i and the localization length £ = Ml, for metallic conduction 12]. The recent 
appearance of accurate persistent current measurements in diffusive electron systems 3], 4 ] 



has been followed by a renewed interest in theoretical studies of such systems 13l. Il5j ]. In 
particular, Ref Jjj] addresses issues such as ensemble averaging and interchannel correlated 
disorder which are related in a broad sense to aspects of this paper. 

The paper is organized as follows. In Sect JIT] we introduce the double-chain Anderson 
model for the disordered ring and derive an exact eigenvalue equation in terms of the transfer 
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matrix in the simple case of pairwise identical random site energies on the two chains. In 
Sect lHII we develop a second order perturbation theory for the two-channel ring eigenvalues. 
In Sect JIV Al we obtain the dominant flux-dependence of the persistent current from the 
eigenvalues averaged over the disorder and in IIVBI we compare our results with the fig. 4 
of {9 1 and discuss the periodicity of the persistent current. Some additional remarks are 
presented in IVl 



II. TWO-CHAIN CYLINDRICAL RING MODEL AND EIGENVALUE EQUA- 
TION 

We consider a vertical cylindrical strip threaded by a flux <p, constituted by two super- 
posed circular tight-binding chains (1 and 2) described by the Anderson model. The system 
obeys the usual tight-binding equations (with a = jrr-) 



-e>* +1 - e- i( Vti = (E- e ln )<p l n - h , (3) 

-e ia <p 2 n+1 - e-**<&_ x = (E- e 2n )<p 2 n - h n = 2, 3, . . . , N - 1 , (4) 

-e ia <pl-e- ia ( p 1 N = (E-e ll ) ( p{-h ( pl , (5) 

-e^(p\-e^ a (p 1 N _ 1 = {E-e 1N )( P 1 N -h^ N , (6) 

e-* a <f N = {E-e 3l )rt-h<p\ , (7) 

-e ia ^\-e-^^ N _ l = {E-e 2N )^ N -h^ N , (8) 

involving one-orbital atomic sites, n = 1,2,... N of spacing a per chainj^, S|. Here 
ip l n denotes the amplitude of an eigenstate wavefunction at site n on chain i (i = 1, 2), £j n 
and E denote the atomic energy of site n on chain i and energy eigenvalues in units of 



minus a constant nearest-neighbour hopping parameter. Finally, h stands for the ratio of 
the hopping parameter between a given site on one chain to its vertical nearest-neighbour 
on the other chain and, minus a constant parameter of nearest-neighbour hopping on the 
individual chains. 

We now transform the Eqs. (J3HHJ) for the two-chain cylindrical strip (wire) by defining 
a basis of independent channels for wave transmission in the absence of disorder. These 
channels are defined by amplitude bases 
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in which the interchain hopping terms in (J2HH]) are diagonal 
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By transforming Eqs (J3HHJ) in the channel bases (fTOl) we obtain successively 



(9) 



(10) 



e^ +1 + e--^-i 



E — h — \{£\ n + £2n) |( £ 2n — Sin) 

2"( £ 2n — £ ln) E + h — \{Si n + £2n) 



E -h - |(en + £ 2 i) 



j( £ 21 — 



5(621 — £ll) + ^ — |(£ll + £ 17v) 

E — h — ^(ein + S2n) \{.£in — £ ijv) 

|(£2JV — £ liv) E + h — ^(Ein + £ 27v) 



(ii) 



(12) 



(13) 



The secular equation for the eigenvalues of the two-chain cylindrical ring are obtained from 
the boundary condition which follows from rewriting Eqs (11 1H13[) in terms of transfer matrices 
P n defined by 



e~ ia P n 



V^-i/ 



n = 2,3,...JV-l , 



(14) 
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yields the secular equation for the energy eigenvalues of the ring: 



(15) 
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4>o 



n=l 



where 



(16) 
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-(E + h) + |(ei„ + e 2 . 



\ 



/ 



(17) 



While the solution of the eigenvalue equation (I16til7j) required for the study of persis- 
tent currents in the cylindrical ring is generally complicated, it simplifies considerably by 
assuming a certain short-range correlation between the site energies on one chain and those 
on the other chain. Thus, in addition to assuming the random site energies E\ m and Ei m 
(rn — 1,2, ... N) to be uncorrelated i.e. 



pm& qn) ^Q^m,n^p,q j (1^) 

as in the usual Anderson model, we impose here that the energy of a site on one chain and 
the energy of the adjacent nearest neighbour on the other chain coincide in any realization 
i.e. 



1,2, 



N 



(19) 



We recall that models with correlated disorder became popular in the context of localization 
in linear chain systems when Phillips and coworkers 16j first showed that they lead to the 
existence of delocalized electron states in one-dimensional disordered systems. In particular, 
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Sedrakyan and Ossipov[17| have recently studied delocalized states in a disordered two-chain 
ladder model with a interchain correlation which is analogous to f fT9|) . For completeness sake 
we mention that the effect of short-range disorder correlation on persistent current has been 
studied recently by us in one dimension, using a simple generalization of the Anderson 
model 18j in which one defines a subset of pairs of nearest-neighbour sites whose energies 
are correlated while the energies of all other pairs are uncorrelated. 

With the correlation (1191) the transfer matrix P n reduces to the simplier block-diagonalized 
form 

P n = P ln ®P 2n , (20) 

where 

f-(E-h)+e n -e' ia \ ~ f-(E + h)+e n -e" ia \ 
Pm= 1 ' ,P 2 n= 1 ' , (20a) 

y e ia J y e ia J 

with e m = ei m = E2m, m = 1,2, . . . N denoting the common value of the random energies 
of adjacent nearest-neighbour sites along the ring. Similarly, under ( !T9|) . the ring transfer 
matrix 

N 



R=l[P n , (21) 



n=l 



with matrix elements Rij, reduces to the block-diagonal form 

R = Ri © R 2 , (22) 

with 

r 1= ( r » r A,% = ( r " r A , (22a) 

\ R21 R22 1 \ R43 R44 1 

where det R\ = det R2 = 1 since R\ and R2 are products of unimodular matrices. 
Using the results (l2"0T 20a). fj2ip and fl22l-22a) the equation ffT6|) factorizes in the form 



l-e 



~(R n + R 2 2) + e ^ detRx 



1 - e~ l ^~ (^33 + R u ) + e~ 1 ^ det R 2 



, (23) 



which reduces to the following final equation for the eigenvalues of the ring in terms of 
partial traces of the transfer matrix 



R u + R22 - 2 cos 



27T0' 



R33 + -R44 — 2 cos 



2vr< 



(24) 



We conclude this section by emphasizing that the ring transfer matrix ( 1221) defined by 
block- diagonal matrices ( 1T7|) (for £i n = e 2 n) constitutes a microscopic realization of the 
transfer matrix of independent conducting channels as envisionned in the phenomenological 
picture of a multichannel wire composed of independent conducting channels. In any case 
the solution of the eigenvalue equation fT24"]) is considerably simplier analytically than the 
solution of ( fTBl for uncorrelated site energies. 



III. WEAK DISORDER PERTURBATION THEORY OF ENERGY EIGENVAL- 
UES 

The energy eigenvalues from which the persistent current may be obtained are related to 
the eigenvalues of the transfer matrix ff22l -22a). In order to study the effect of a weak disorder 
by a perturbation theory we first determine the eigenvalues in the absence of disorder. In 
this case 



R = P N , with P = Q 1 ®Q, 



where 



Qi 



-(E-h) -e 




(25) 



(25a) 



In order to obtain the elements of P we diagonalize P by means of a similarity transfor- 
mation. The eigenvalues of P are 



ai 



A: 



—(E -h)± ^{E-h) 2 -A 
-(E + h)±y/(E + h) 2 -4: 



(26a) 
(26b) 



which correspond to unperturbed energy bands in terms of wavenumber variables s\,S2 
defined by 
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Xi = e ±4Sl (27a) 
Xl = e ±lS2 , (27b) 



with 



E = E h ( Sl ) = h-2coss 1 , (28a) 
E = E_ h {s 2 ) = -h-2coss 2 (28b) 

where we assume h > 1 and refer to the si- and S2-band as the upper and the lower band, 
respectively. The similarity transformation which diagonalizes P in terms of the eigenvalues 
( 1271) is defined by the eigenvector matrix 



U = U 1 ®U 2 (29) 

with 

( e i(si-a) e ~i{si+a) \ 
^ ^ \,i = l,2 . (29a) 

Using the relation P m = U(U- 1 PU) m U- 1 we finally obtain, from fl27ati27bl) and (|29l29a), 

P m = Q™ © Q™ (30) 

with 

~ 1 / sinfm + l)s{ —smmsie~ ia \ 

QT = —\ \,i = l,2 . (30a) 

smsj l sinmsje"* — sin(m — l)sj / 

In the absence of disorder we insert ff25l) . fl30|) and (30a) in ff24|) and obtain for the wavenum- 
bers in the energy bands fl28aH28b|) 



8l = 8l (k) = ^(k + ^,k = l,2,...,N , 
S2 = S2 (k') = ^(k , + ^j,k , = l,2,...,N 



(31a) 
(31b) 
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Next we develop a perturbation theory to discuss the effect of a weak site-energy disorder 
on the eigenvalues of the two-channel ring threaded by a flux cf). This generalizes our earlier 
treatment js ] of the effect of disorder on persistent current in a one- dimensional ring. The 
latter |8j has also been applied recently for studying the effect of short-range correlated 
disorder on persistent current in one dimension 18]. 

We now wish to solve fl24"]) for weak random fluctuations of the site energies. To this end we 
expand the transfer matrix to successive orders in the £j in f lT9|) and obtain the corresponding 
first and second order corrections for the perturbed energy levels of the (coupled) two-chain 
system of the form 



k,k' = 0,±1,±2 



(32) 



where 



E i% (k, k') = -2(cos Sl + cos 8 2 ) • (32a) 
The transfer matrices P n in ffTTj) are split into unperturbed and perturbed parts 



P n = P + V n = P ~ 



AE - -(e ln + e 2r . 









e © e, e = 













where 



AE = + E {2) + 



and the matrix product Yln=i ^ n ^ s expanded to quadratic order in the quantities 



[£ln + ^2n) 



This yields 



is 



(33) 



(34) 



AE- 



N 



N 



R=l[P n = P N +J2 P m ~ l V m P 

n=l m=l 

N n-1 

p m ' 1 v m p n ~ m ' 1 v n p N - n + 



N-m 



(35) 



n=2 m=l 



Note that the expansion in powers of V n in (|35|) goes beyond a systematic expansion to 
successive orders in the site energy perturbation since e.g. the terms linear in V n , besides 
being linear in e n , depend on all orders of the energy level perturbation via AE. Nevertheless 
( 1351) leads to a systematic determination of the successive order perturbations of the energy 
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levels since the solution of the eigenvalue equation (|24|) imposes that the disorder effects in 
the partial traces of the transfer matrix must vanish. The evaluation of the first and second 
order corrections in the energy levels (1321) proceeds by performing the matrix products 
entering in fl33|) . in order to obtain the diagonal elements of the ring transfer matrix in the 
partial traces in (EMI) . 

The contribution of the linear correction term in (155)) to the trace of R is given by 



TyR 



1 1 
+ — 



N r 1 



27T0 



sin si sin s 2 J (po *-^f. ' ' 2 



+ 0(e< 



(36) 



using (150)) . (30a) and (133)) . Since according to (JMJ) eigenvalues of the disordered two-channel 
ring correspond to values of the partial traces which are independent of the disorder, it follows 
from ( 136)) that the first order effect in the eigenvalues in the correlated site energies model 
is 



E sll 2 = ^y5Z^ lm + £ ^ 



(37) 



except for flux values equal to integer multiples of 4>q/2. Note that these flux values are 
precisely those at which the free particle spectrum of the ring is degenerate [7]. 

At quadratic order in the site energies, on the other hand, there are two types of contri- 
butions in Eq. (155)) . R = Rl + R", the first one 



N 



N—m 



(38) 



m=l 



and the second one 



N n-l 

n " = E E^£« - " £ rn)P m -\e © e)P^ m -\e © e)P N ~ n . (39) 

n=2 m=l 

The explicit evaluation TrR', using (150)) and (30a), yields 



Tr^' = E?l 



sin 2 Si 



A^cossi cos A^si 



3 sin A^si 

4 sin 3 si 



same with si — >■ s 2 



(40) 



Similarly, the evaluation of Tri?" leads to 
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n=2 m=l 



sin(n — m)sism(N + m — n)si 

2 hsame with si -> s 2 

sin Si 

(41) 

The condition TrR' + TtR" = which follows from ( l24"j) thus yields the final result 



1 /\r »r 3siniVsi\ . , 1 1 

Jy cos si cos TV Si — T^r^i — + same with s\ — > s 2 



sin Si \ 4 sin s\ 

sin(n — m)si sin(iV + m — n)s\ 



N n-1 

sm (/; — ;/)). si sm ( \ ■■; in — n)s> 

+ same with si — > s 2 



n=2 m=l 



sin 2 si 



(42) 



which will be exploited in the following section for discussing the effect of the disorder on 
the persistent current. 

IV. EFFECT OF DISORDER ON THE PERSISTENT CURRENT 
A. Calculation of the persistent current 

In this subsection, like in our previous work we discuss the explicit form of the 

persistent current obtained simply from the energy levels of the ring averaged over the 
disorder. Indeed, the more accurate approach which consists in finding the persistent current 
(using (2a) from Eq. ( j4"2l before averaging over the disorder leads to results whose further 
analysis is intractable analytically. By averaging (142p . over the disorder, using ( JT8lfT9l) . we 
get, after some simple algebra, 

(7V-l)cos2iVsi sin(7V-l)si i ■, , 

„ I 1 7-3 - -3 + same with si — > s 2 

2 \ sirr si sin' 5 si 1 ' 



N cos ai cos Ns, _ SsmNsi ] + game with ^ _^ ^ 



{Eifjk,k')) = - £ ^ 1 — L , (43) 

us si cos Ns\ 3 sin Nsi 

sm A s\ 4 sin 5 s\ 

which reduces, for large N, to the form 

(2) 4( ^l^f + same with si ^ s 2 

( E s in ) = "J 1 J osZfNs, ~ same with ~ ~ ^ ) , ( 43a ) 
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up to corrections of order 1/JV. Using f)31aH31b]) this expression may be rewritten exactly in 
terms of an important fixed flux-dependent prefactor which is common to all energy levels 
(E slS2 (k, k')) of the ring: 



(*%,(*. = C 4 Pet - * co. ...... : ^ - + : • («) 



4 V cos ^ </>o 7 cos Si sin 2 s 2 + cos s 2 sin 2 si 

<Po 

Note that this expression is not valid for flux values equal to odd integer multiples of 0o/4 
at which the perturbation theory clearly breaks down. 

The effect of the disorder on the persistent current of an electron in an energy level f )32|) of 
the ring defined approximately from the averaged second order energy perturbation, namely 

(2) d(Egl 2 ) (A) 



is given by 



e(M'»^-e^(sin^)(2 + - 1 ¥ )/. 1 « +0 (i) , (46) 
with 

f tu u>\ sin 2 si + sin 2 s 2 . . 

fs lS2 {k,k) = — ■ r^— , (46a) 

cos Si sin s 2 + cos s 2 sin S\ 

whose differentiation in ( I45p yields the terms of order 1/N which have been ignored in 
(f46|) . Clearly, for 0/0o N/2n the coefficient of f SlS2 in (T46l) represents the dominant flux 
dependence of the disorder contribution, 1^ = ])2 k k / Isfs 2 (k,k',(j)) in the total persistent 



current in the two-channel system 191 ] . In particular, this coefficient reveals a strong decrease 
of the persistent current due to the disorder for fluxes approaching the value 0o/4 from either 
side (the value 0o/4 itself lying obviously outside the domain of validity of the perturbation 
theory). The explicit (numerical) evaluation of I^ 2 \ say in the typical case of a half-filled 
upper energy-band S\ and a completely filled lower band s 2 lies beyond the framework of 
our analytic treatment and will not be pursued further. 



B. Persistent current periodicity, single- and two-channel effects in simulation 
results for a multi-channel ringQ] 

There exist two fundamental periods in the persistent current problem, namely the 
nearest-neighbour hopping period N<j) in the hamiltonian of Eqs. (J3HHJ) and the revolution 
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period 0o for an electron on the ring subjected to the flux 0, which affects the boundary 
condition ([Q). The periodicity iV0o is responsible for the period 0o/2 of the free particle 
persistent current averaged over even and odd numbers of electrons in a one-dimensional 



7. 



ring [3, |8| as recalled below. On the other hand, the periodicity 0o/2 is also obtained exactly 
for the effect of the disorder in a one-dimensional ring. Since, however, the 0o/2 period- 
icity of the averaged free particle persistent current was not explicitely mentioned in 
we briefly discuss it here, along with the influence of the disorder in one dimension. These 
results are important for our qualitative interpretation below of salient features of numerical 
simulation results of Bouchiat and Montambaux (fig. 4 of [9[) for the effect of weak disorder 
on persistent current in a ten-channel (cylindrical) ring, in terms of single- and two-channel 
contributions obtained in Ref. [8| and in the present paper, Eqs ( l4"6"l 46a). 
Let 

i° = (iL) = l(iLn + C?) , (47) 

be the persistent current in the unperturbed half- filled upper tight-binding band in (32a), 
averaged over even and odd numbers of electrons iVe|lO|. The currents 7° ven and J° dd are 
given by 0, Q 



lLn--Io -r- 1 . ( 48a ) 



^0 

/ O ° dd --/o^,0<^<i , (48b) 

00 0o 2 

and are illustrated by the figures 3(a) and 3(b) of [7j (whose legends should be inverted 



181]) . Their determination by summing contributions from the occupied energy levels below 
the Fermi level is easily visualized from the energy level spectrum as a function of flux shown 
in fig. 2 of Q. The current in ( 1471) is thus 

whose periodicity is Here Iq = ^ is the persistent current of an electron at the Fermi 
level and Vp = j[-jr is the Fermi velocity. Next we recall the effect of a weak disorder on 
the magnitude of the persistent current in a one- dimensional ring given by Eqs (I49H50P of 



After averaging the effects obtained in [8] for even and for odd numbers of electrons as 
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in (149|) . we obtain the total persistent current in the disordered half- filled one- dimensional 
ring: 



/=-/n 



f - 



Nel 



Isin 2 (27r0/0 O ) 



,0< 



t 1 
60 2 



(50) 



which also has periodicity 0o/2. We note, incidentally, that ( 1501 provides a simple qualita- 
tive interpretation of the numerical results of Trivedi and Browne ll|] for persistent current 
in one-dimensional disordered rings averaged over even and odd numbers of electrons (fig. 



4 of 



Uj and fig. 8 of (2]). 



At present we wish to compare our results for persistent current averaged over electron 
numbers in the single- and two-channel rings with dominant features shown by the flux- 
dependence of similarly averaged persistent current in a ten-channel ring in fig. 4 of Bouchiat 
and Montambaux {J. The 10-channel disordered ring is described by the Anderson-model 
using a rectangular distribution of the site energies of width W = 0.2 (in units of the nearest- 
neighbour hopping parameter) [20] and the persistent current is averaged over numbers of 
electrons 150 < N e < 250 in the range < <p < </>o/2. 

The pertinence of a comparison of our results for one- and two-channel systems with 
the multichannels results of fig. 4 of [q| is due to the fact that the latter were found to 
be independent of the number of channels^]. On the other hand, the inclusion of results 
for the two-channel system besides those for a one-dimensional (single-channel) ring for the 
comparison with the fig. 4 of 9J enriches our model through the novel effects related to the 
interchain hopping. 

We now come to the detailed comparison of our perturbation results for persistent currents 
in one- and two-channel disordered rings with the numerical results for a ten-channel system 
in fig. 4 of Ref. First we consider the free particle persistent current ( 1491) which shows, 
in particular, the period halving effect resulting from averaging the current over even and 
odd numbers of electrons in fig. 3 of [4]. This expression provides clearly a good "zeroth 
order" approximation of the overall current profile in fig. 4 in [9j. Next, the disorder effect 
in the persistent current (j50p averaged over numbers of electrons, for a one- dimensional ring 
accounts qualitatively for the relatively large current dip near = and a corresponding 
peak near <fi = <po/2 due to disorder, in fig. 4 of 9|. Finally, the disorder contribution 
j(2) _ k , Isil 2 in the persistent current in the two-channel ring involves the dominant 
fixed flux- dependent overall factor in (|4"6j) . which leads to a strong decrease with respect to 
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the uniform decay ( 149]) for — > 0o/4 from below and to a similar strong decrease with respect 

to the uniform increase of current for — > 0o/4 from above. These current reductions are 

in line with a corresponding strong reduction of current with respect to a uniform decrease 

in fig. 4 of 9] (which passes through zero at = o /4) for — > 0o/4 from below, starting 

from = 0, 216 O up to = 0, 233 O , and with a corresponding reduction of current 

fl 

with respect to a uniform increase of unperturbed current in fig. 4 of [9j] for — > 0o/4 
from above, starting from the value ~ 0,275 0o down to ~ 0,254 0o. Finally, in our 
reference to the fig. 4 of [9J we have excluded the domain 0, 233 < < 0, 254 0o from our 
discussion not just for convenience but mainly because the perturbation theory breaks down 
at = 0.25 0o- On the other hand, our identification of the most significant features of 
the persistent current in multi-channel system with single- channel and coupled two-channel 
disorder effects suggests the possibility of identifying further finer features of the fig. 4 of 
Bouchiat and Montambaux 9J with 3-4, ...-channel effects. This might suggest an analogy 
between the persistent current problem in multi-channel systems and e.g. the interactions in 
a typical many-body system such as a molecular crystal in which the many-body interaction 
between the molecular species may generally be decomposed into convergent pair-, triplet-, 
... interactions. 

Finally, we turn to the discussion of the periodicity of the persistent current. On the basis 
of their extensive numerical simulations for a weakly disordered 10-channels ring Bouchiat 
and Montambaux 9j have found that the persistent current averaged over a range of electron 
numbers has periodicity 0o/2. However, the physical origin of this period halving remained 



unclear 2l|. As shown by ( 1501) the analysis of 8( leads to period halving of the persistent 
current in a disordered one-dimensional ring while the discussion below for the two-channel 
ring indicates that besides the periodicity 0o/2 there exists another approximate periodicity, 
0o, which is induced by the disorder. 

As mentioned earlier we consider the typical case of a two-channel ring where the lower 



energy-band ( I28bl I31bl) is completely filled and the upper band (I28all31al) is half-filled, with 



N e = N/2 electrons. The filled lower band carries no free particle current (vanishing fermi 
velocity) and the total current in the upper band is approximately given by 

(26 1\ elite 2tt0 / 1 \ 1 

I — — Iq [ —-— — — I 0— — sin — -— ( 2 + jo<-7-<- . 51 

°V0o 2 J 2h O V cos 2 ^;' 0o 2 K J 



16 



The first term is the free particle current ( 149]) averaged over even and odd numbers of 
electrons in the upper band. The second term determines the dominant decrease of the 
persistent current 1^ = ^2 k k , li\l 2 (k } k' } 0) in the two-channel ring which is induced by the 
disorder. It is given by by ignoring the weak perturbation arising from summation 
of f Sl s 2 (k,k') over the occupied individual energy levels averaged over electron numbers, 
whose magnitude is typically of order 1 at ~ 0o/4, relative to the larger effect of the 
dominant flux dependence in (15T1) around this flux value. This result for the order of mag- 
nitude of J2 k k , f SlS2 (k, k') is inferred from the magnitude of the corresponding simpler sum, 

Ylik fs{k) = Yl V s i ns ) s — + -^j ,k = 1, 2, ... N, which arises in the study of the 

effect of disorder on persistent current in a one-dimensional ring [8| and whose approximate 
evaluation yields results linear in 0/0o both for even and for odd numbers of electrons, for 
< (p/4>o < 1/2. Finally we note that the above estimate of the persistent current near 
— 0o/4, which focuses on finding its order of magnitude, does not necessarily yield the 
correct sign of the current. 

The existence of two distinct components of persistent current, a free par ticle component 
of periodicity 0o/2 and a component of periodicity 0o due to the disorder [22| may reflect the 
presence of two distinct periods in the tight-binding equations (J3}|8j) (period N<j) ) and in 
the boundary condition (JT]) (period 0o)- Note that in the single-channel (one-dimensional) 
case the disorder effect itself in the persistent current before averaging over particle numbers 
involves the periodicity 0o/2 (see Eqs. (49-50) of 



V. CONCLUDING REMARKS 



We conclude with a general remark concerning the possibility of interpreting further 
specific features of the numerical results for persistent current averaged over numbers of 
electrons in many-channel disordered systems in fig. 4 of Bouchiat and Montambauxjol in 
terms of disorder effects in few-channel systems, say M = 3 or 4. 

The definition of independent channels in a non-disordered M-chain tight-binding cylin- 
der threaded by a flux involves in a first step the diagonalization of the interchain hopping, 
t, leading to the energy eigenvalues [23] 
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M+l 

for the channels n = 1, 2, . . . M, where & n = 0, ±1, ±2 . . .. For M = 2 this expression reduces 
clearly to ( 1281) with the definition (13T1) of wavenumbers. However, for larger M the discussion 
of the effect of disorder becomes rapidly very complicated analytically. Nevertheless, specific 
features in the persistent current, such as the existence of weaker peaks and/or dips at special 
flux values, which are related to the interchain hopping may be expected in analogy to what 
we have found for the single-chain- and for the two-chain cylindrical ring model studied 
above. 
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